Small Fermi energy and phonon anharmonicity in MgB 2 and related compounds 
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The remarkable anharmonicity of the E^g phonon in MgB2 has been suggested in literature to play 
a primary role in its superconducting pairing. We investigate, by means of LDA calculations, the 
microscopic origin of such an anharmonicity in MgB2, AIB2, and in heavily hole-doped graphite. We 
find that the anharmonic character of the E2 g phonon is essentially driven by the small Fermi energy 
of the a holes. We present a simple analytic model which allows us to understand in microscopic 
terms the role of the small Fermi energy and of the electronic structure. The relation between 
anharmonicity and nonadiabaticity is pointed out and discussed in relation to various materials. 

PACS numbers: 74.70.Ad, 63.20.Ry, 63.20.Kr 



The report of superconductivity at T c — 39 K in MgB2 
jl| has raised great expectations about metal diborides 
MB2. There is indeed no reason to believe that MgB2 
represents the highest-T c compound within this family. 
Since the very beginning a generic consensus about the 
electron-phonon (e-ph) nature of the superconducting 
pairing has prevailed [g, |||, although some purely elec- 
tronic models have been proposed Q. The precise ori- 
gin of such a high- T c superconducting phase is still un- 
known. Recently the strong anharmonic character of the 
in-plane Ei g phonon mode and its possible correlation 
with the high T c value in MgB2 have attracted a consid- 
erable interest |5|, fj], || [|. Here we attempt a sim- 
ple theory of such a strong anharmonicity and test its 
predictive power on related compounds. With this per- 
spective we have performed first-principles calculations 
of the band structure and lattice properties of MgB2, 
A1B 2 , and of a hypothetical hole-doped graphite. We 
identify the small value of the Fermi energy for the holes 
in the a band, with entire portions of the Fermi surface 
disappearing upon Ei g distortion as the funda- 

mental origin of anharmonicity; a simple way of modeling 
the effect of distortion on the band structure confirms our 
finding. It has already been pointed out in literature that 
MgB2 resembles in many ways graphite 



if VA tta 



From the structural point of view MgB2 is formed by 
graphenc like layers of B spaced by planes of Mg atoms. 
The point group symmetry of in-plane boron phonon 
modes of MgB 2 and in-plane phonons of graphite is thus 
the same fl, and the difference in frequency is related 
to the different strength of B-B and C-C bonds. Out of 
the whole vibrational spectrum, a large interest has con- 
verged towards the Ei g mode in MgB2, which involves 
only in-plane boron displacements. This mode has been 
shown to have an extremely strong coupling with the in- 
plane a bands §, |, |, |, 0, [§, § , which in MgB 2 provide 
conduction holes; a relation to the high superconduct- 
ing temperature T c was naturally suggested. This idea 



was also supported by the negligible partial isotope coef- 
ficient on T c associated with the Mg atomic mass . In 
this context the strong anharmonicity, a unique property 
of the E-^a phonon within the MgB2 vibrational spec- 
trum (HI £j, acquires an obvious importance. The elec- 
tronic structure of MgB2 shares strong similarities with 
graphite. In both materials one can identify strongly 
two-dimensional a bands, almost entirely derived from 
B (C) s and p x ^ v orbitals, plus bonding and antibonding 
7r bands with three-dimensional dispersion and mainly 
Bp- Mg (C Pr ) character. The most important difference 
between MgB2 and graphite is the position of the Fermi 
level In undoped graphite n cuts the tt band struc- 
ture just in the middle, about 3 eV above e^ op , the top of 
the bonding a bands, which are thus completely full Jl5| . 
In MgB2 the combined effect of additional magnesium 
layers, different ionic charge between boron and carbon, 
and valence-charge transfer from magnesium to boron 
layers, yields a different arrangement of a and 7r bands. 
The resulting Fermi level /1 still cuts the 7r band struc- 
ture somewhere, but is now ~0.5 eV below the top of the 
a bands, which therefore, in MgB2, give a sizable hole 
contribution to the Fermi surface jllL [13| , not present in 
graphite. This important difference was almost immedi- 
ately pointed out more recently the large splitting 
and shifts which the a bands undergo upon typical E2 g 
phonon displacements were suggested as the likely source 
of anharmonicity for the E2 g phonon in MgB2 ||. We 
are going to make this statement more precise. We claim 
that neither the presence at the Fermi level of the a bands 
nor their strong coupling to the Ei g phonon are suffi- 
cient to induce anharmonic effects: it's the small Fermi 
energy associated, in the unperturbed crystal, to the cr 
conduction holes (e^ op — fj, ~ 0.45 eV) which makes the 
difference. On these grounds we may expect anharmonic 
effects to be strong in other materials with small Fermi 
energies (and sufficient e-ph coupling), and also conjec- 
ture a relation between anharmonic and nonadiabatic ef- 
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fects. We present local Density Functional Eg] calcula- 
tions, based on Martins- Troullier pseudopotentials fTij] 
and the ABINIT code |l8|], which support this picture; 
the resulting bands and E 2g frozen-phonon energies are 
well understood in terms of a simple model, discussed in 
the second part of this letter. We studied the E 2g phonon 
for MgB2, AIB2, graphite and a hypothetical hole doped 
graphite where one electron is missing from each carbon 
atom. Experimental lattice parameters were used as an 
input for MgB2, AIB2 and graphite, while for hole-doped 
graphite the lattice parameters were left at the experi- 
mental value of graphite. For the purpose of this study 
the use of GGA's and an extremely accurate a priori de- 
termination of equilibrium lattice parameters are not an 
issue. Great care in the k-space integration is, instead, an 
issue 0, since, in some cases, entire portions of the Fermi 
surface disappear upon distortion; we use two 15x15x10 
shifted Monkhorst-Pack grids in the Brillouin zone |L8| . 
We should also specify that, between the two possible 
eigenvectors for the E 2g phonon, we only show results 
for the one labeled E 2g (b) in Ref.||, whose energy is, by 
symmetry, an even function of the displacement; phonon 
displacements up to 0.05 ~ 0.1 A were considered. The 
first important observation is that the E 2g anharmonic- 
ity is completely absent in AIB2 , whose a bands undergo 
equally large splittings and shifts as MgB2. As shown 
in Fig. |l|, our frozen-i^g-phonon calculations reproduce 
the large anharmonicity found for MgB2 (black squares) 
H |(| , but predict no such effect for A1B 2 (empty squares) , 
whose energy remains proportional to the square of the 
phonon displacement for all the displacements under con- 
sideration. Our finding is consistent with the experimen- 
tal fact that the E 2g phonon line is very broad in MgB2 
but not in A1B 2 §). 
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FIG. 1: Energy AE associated to an E^ g phonon displace- 
ment of amplitude u, plotted as a function of u 2 . For each 
material this energy is divided by 02 (in the inset, units of 
eV/A 2 ), the quadratic coefficient of a polynomial best fit 
AE ~ a,2U 2 + a,4u 4 + .... On both axes the units are thus 
A 2 , and harmonic phonons collapse on a single straight line 
y = x. The solid lines result from our model (Eqs. ^ and 
the corresponding parameters are shown in Table [l| 
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FIG. 2: Upper panels: MgB2 electronic bands without (left) 
and with (right) an i?2 9 phonon distortion of amplitude u = 
0.05 A. Lower panels: same as upper panels, but for AIB2 
and u — 0.05x(o,aib 2 / a MgB 2 ) A. The a bands are marked as 
thicker lines. Their splitting upon E^g distortion is equally 
large in MgE>2 and AIB2, but in the latter (lower panels) they 
are always below the Fermi level (/i = in all panels). 



We suggest that the different behavior of AIB2 be sim- 
ply related to the fact that, both before and after the E 2g 
distortion, its a bands are near, but completely below, 
the Fermi level /1; unlike MgB2, their electronic occupa- 
tion remains unchanged upon distortion (see Fig. |^) . In 
MgB2, instead, the top of the a bands is above the Fermi 



3 



energy but, upon distortion, the lower splitoff band com- 
pletely sinks below it, thus changing its occupation. Be- 
sides the amount of the shifts and splittings, the exact 
position of the a bands before and after the E 2g dis- 
tortion is thus a crucial ingredient for its anharmonic- 
ity. This is confirmed by artificially moving the top of 
these bands, e^ op , w.r.t. the Fermi level /i in graphite. In 
true graphite the a bands also undergo large splittings, 
but they are already well below \i both before and after 
the E 2g distortion; here our frozen-i^g-phonon calcula- 
tions find no anharmonicity (Fig. [IJ empty dots; "gr." 
stays for graphite). But if, by adding a uniform neutral- 
izing background fll3| , we remove one electron per car- 
bon atom (two electrons per cell, black dots in Fig. [I]) 
from graphite, thus shifting e^ op , the top of its a bands, 
back to the "optimal" position (~ 1 eV above /i in the 
undistorted crystal), then, upon typical E 2g distortions, 
the lower splitoff a band sinks below /i, and we recover 
strong anharmonic effects, shown in Fig. [I] (black dots). 
To clarify the origin of this anharmonic behavior we have 
traced back the effect of the E 2g frozen-phonon distor- 
tion to the electronic structure. For all materials we find 
that the main effect of the lattice displacement is a lin- 
ear energy splitting of the a bands. The it bands are, 
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TABLE I: Three inputs for our total energy model, Eqs. (y) 
and (^), extracted from our LDA outputs. The remaining 
two parameters No- and N„-, needed only when e CT ° p > /x, were 
adjusted to yield the best fits shown in Fig. [jj Their optimal 
values (N CT =0.11, = 0.39 for MgB 2 ; N CT =0.07, N^ = 0.30 
for graphite+- 1-) fall in a physically reasonable range, in spite 
of our oversimplified density of states (see text). 

instead, only weakly modified. The effects of lattice dis- 
tortion on the band structure, in a relevant energy range 
around the Fermi level [a, can be thus schematized, to a 
good approximation, as <5e CT (k) ~ ±g|u|, with opposite 
signs for the two different a bands. The value of g de- 
duced from the LDA results differs significantly between 
the borides (MgB2, AIB2) and graphite (with or without 
doping), while it is almost constant within each class of 
compounds (see Table [j]). ^From the comparison of the 
band structures in Fig 70 we can identify two represen- 
tative cases. For AIB2 (bottom) the bonding a bands, 
at zero distortion (left), are completely below the Fermi 
level [i, and remain there upon distortion (right); their 
energy splitting, induced by the E 2g lattice displacement, 
does not change either their occupation or the topology 
of the Fermi surface, which never acquires a a sheet; no 
a band crosses /i at any displacement, and the Fermi 
surface is exclusively dictated by the 7r bands, which, 
compared to the a bands, undergo only minor changes 



upon displacement. Besides AIB2, this is also the case of 
undoped graphite (not shown) . An entirely different situ- 
ation is found when, on one hand, the top of the bonding 
a bands at zero distortion, e^ op , is above the Fermi level 
/j, (so that in the perfect crystal the Fermi surface has a- 
hole-like cylindrical sheets p]), but, on the other, the 
energy splitting of the a bands is large enough to drive 
one of them completely below /j, upon distortion. This 
is the case of MgB2 (top Fig. ||), and also of heavily 
hole-doped graphite (not shown). In both cases, at some 
critical phonon displacement, the number of Fermi sur- 
faces associated to the a bands changes (one of the two 
cylindrical sheets around the T-A line disappears M); 
beyond that point, since the total number of electrons is 
conserved, larger displacements will imply a qualitatively 
different behavior, due to the reshovcling of electrons be- 
tween a and tt states. To gain further insight, we present 
a simple model of the E 2g anharmonicity which seems to 
represent well both types of situations. We first consider 
the system with no distorsion (u = 0). The electronic 
band energy can be written as: 

%=0)=2^e i (kK(k) + 2^ £ ,(kK(k), (1) 

k,i k 

where e^(k) represents the dispersion relation of the two 
a bands, e^k) takes into account the remaining ir bands, 
and rii(k), n^Ck) are the corresponding occupations (the 
factor 2 is for the spin degeneracy). The LDA bands 
teach us that the most important effect of an E 2g distor- 
tion u ^ is an almost linear splitting (i.e. proportional 
to u) of the two a bands around the T-A line. This can 
be roughly modeled by a linear e-ph Jahn-Teller-like 
coupling of the E 2g phonon to the a bands; the small 
coupling to the ir bands is neglected altogether. The re- 
sulting total energy at u^0 is: 

E(u)= 2 5^ei(k)ni(k,u) +2^ f7r (kK(k,ti) 

k,i k 

+ 2gu}2[MKu)-n 1 (k,u)} + ^^u 2 (2) 

k 

In Eq. (||) the electronic band energy, in the presence of 
a phonon displacement u, was split into the sum of three 
terms: unperturbed a and ir bands with it-dependent 
occupation (first two terms in the r.h.s), plus a linear e- 
ph coupling of the E 2g mode with the a bands. The 
last term is an effective elastic energy. Far from the 
Fermi level our bands are not realistic (see below), so we 
lump into this term both the bare ion-ion repulsion and 
those electronic effects which are missing from our model 
bands. The occupation number can be self-consistently 
calculated: ni(k, u) — /[ei(k) — gu — fi(u)], ri2(k,ii) = 
/[e 2 (k) + gu - n{u)}, n ff (k, u) = /[e^k) - fi(u)], where 
/x(it) is the Fermi level in the presence of the frozen 
phonon and f[x] — 9{—x) is the T = Fermi function. 
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Let us consider the representative case of MgB2. For 
the sake of simplicity we assume two parabolic a bands, 
perfectly two-dimensional and degenerate at u = 
[ei(k) = e 2 (k) = e CT (k)]. The corresponding density of 
states (DOS) of each a band will be therefore constant 
up the top of the band e^: N a (e) = N a [e < 4° p ]. From 
now on we conveniently set n(u = 0) = 0, so that e^ op 
now equals e^ op — fi(u = 0), the Fermi energy of the a 
holes in the absence of lattice distorsion. In addition, 
in the energy range we are interested of, we can assume 
iV-Tr , the density of states of the it band, to be just con- 
stant. The Fermi level /i(u) and the total energy E(u) 
in the presence of the frozen phonon distorsion u can 
now be easily computed. The system shows a qualita- 
tively different behavior for two interesting regimes: (i) 
g\u\ < e^ op and (ii) g\u\ > e^ op . Within the assump- 
tions of our model (rectangular DOS), in the regime (i) 
the Fermi level is unaffected by the frozen phonon dis- 
torsion, /j,(u)=fi(u = 0) = 0: the depletion of the a band, 
raised by the distortion, is compensated by an equivalent 
filling of the other a band, lowered by the same amount; 
the 7r bands, modeled by their constant DOS N w , play 
no role. We obtain for the total energy: 



E(u) = E(0) 



^ku 2 -2N„ g 2 u 2 



9\u 



<4° P - (3) 



Eq. (|^) represents a phonon frequency renormalization 
due to the response of the a electrons, E(u) = E(0) + 



Mnj g u 2 /2, with n 



J 2g 



AN a g 2 /M. The harmonic 



character of the Ei g phonon mode is however unaffected. 
Things change in the regime (ii) g\u\ > e^ op . When the 
energy splitting is larger than the zero-distortion Fermi 
energy of the a holes, g\u\ < e^ op , the lower a band is 
completely shifted below the Fermi level. Now this band 
is full, and cannot further compensate the loss of elec- 
trons from the upper a band. Then the only way to 
conserve their total number is to add more electrons to 
the 7r bands, which thus come into play. To obtain this, 
the Fermi level needs to shift, and the dependence on u 
of the total energy, is, in turn, deeply modified: 



H(u) 



top 



) 9\u\ > e 



lop 



E(u) 



= E(0) + -^-u 2 - 2N a g 2 u 2 



(4) 



(5) 



JV ff (2iV g +JV ff ) 



,top\2 



) 2 5H>4 op - 



In our simple model the transition between harmonic and 
anharmonic regime occurs when one band is completely 
shifted below the Fermi level, and does not manifest itself 
as a simple additional quartic term u A . Rather, an over- 
all anharmonic potential results from a simple harmonic 
term up to g\u\ < e^ op which, for g\u\ > e^ op , smoothly 
connects to a shifted parabola with different curvature. 



Such a non-analytic behaviour has to do with the ex- 
treme simplifications of our model, in particular with the 
perfectly 2D parabolic character of the a bands (step- 
like density of states), and on the assumption of per- 
fect degeneracy of the a bands at u — [ei (k) = e 2 (k)] . 
We have checked that, with slightly more realistic mod- 
els, the sharp transition of Eq. Q becomes considerably 
smoother. However, Eq. (|5|) is particularly appealing 
just because of its simplicity, since it depends only on a 
few parameters which can be extracted from LDA cal- 
culations (see Table [jj), thus providing a direct test of 
the model. The results are shown as solid lines in Fig. |l| 
for the compounds considered here. The agreement with 
LDA first-principles calculations is quite good, consider- 
ing the extreme simplifications of our model. In conclu- 
sion, in this paper we have investigated, by means of nu- 
merical and analytical techniques, the microscopic nature 
of the anharmonicity of the E^ g phonon mode in MgB2. 
The results presented here provide a clear evidence that 
the anharmonicity of the E^ g phonon mode in MgB2 is 
induced by the extremely small value of the cr-hole Fermi 
energy. Along this view we can predict a strong anhar- 
monicity in heavily hole-doped graphite. We have shown 
that the anharmonicity of thei?2g mode, which is strongly 
coupled to the a bands, can be considered a signature 
of small Fermi energy; this points out in a natural way 
towards the possibility of nonadiabatic effects. A quan- 
titative description of this situation involves, however, 
quantum-many-body effects (nonadiabatic renormaliza- 
tion of the phonon frequencies Jf9[ ) . While this task is 
beyond the aim of the present paper, different theoretical 
studies already suggest that nonadiabatic effects could be 
responsible for the high T c in MgB 2 |2(], |^|. Note that 
this is entirely different from the initial claim that 
anharmonicity affects superconductivity via the nonlin- 
ear coupling. From this point of view the situation is, 
instead, similar to fullerenes |2^] , which have been re- 
cently shown to reach critical temperatures as high as 
T c = 117K in FET doped compounds @. In this re- 
spect our work suggests new perspectives in the search 
for high-T c materials. In particular, heavily hole-doped 
graphite, which we predict to have small Fermi energy 
and anharmonic E2 g phonon, would be a potential candi- 
date for high- Tc superconductivity. The recent claim of 
T c = 35K in amorphous graphite-sulfur composite sam- 
ples could be related to this scenario p4|. 
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